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We study the effect of insulating oxide substrates on the energy band structure of monolayer 
and bilayer graphene using a first principles density functional based electronic structure method 
and a local exchange correlation approximation. We consider two crystalline substrates, Si02 (or 
Q-quartz) and AI2O3 (a-alumina or sapphire), each with two surface terminations. We focus on 
the role of substrate surface dangling states and their passivation in perturbing the linear energy 
spectrum of graphene. On non-passivated surface terminations, with the relaxation of top surface 
layers, only Si-terminated quartz retains the linear band structure of graphene due to relatively 
large equilibrium separation from the graphene layer whereas the other three surface terminations 
considerably distort it. However, without relaxations of the top surface layer atoms, linear bands 
appear in the electronic spectrum but with the Dirac point shifted away from the Fermi level. 
Interestingly, with a second carbon layer on non-passivated oxygen terminated Quartz, with top 
surface layers relaxation, graphene features appear in the spectrum but sapphire with both surface 
terminations shows perturbed features even with two carbon layers. By passivating the surface 
dangling states with hydrogen atoms and without top layer atomic relaxations, the electron-hole 
symmetry occurs exactly at the Fermi level. This suggests that surface dangling states play a less 
important role than the atomic relaxations of the top surface layers in distorting the linear spectrum. 
In all cases we find that the first layer of graphene forms ripples, much like in suspended graphene, 
but the strength of rippling is found to be weaker probably due to the presence of the substrate. 
We discuss the energetics of both passivated and non-passivated surface terminations with a top 
graphene layer. 

PACS numbers: 71. 15. Mb, 71.20-b, 73.20.-r 



I. INTRODUCTION 

The excitement in graphene research, due to the real- 
ization of table-top high-energy physics experiments and 
the promise to replace silicon in future semiconductor 
chips, is certainly overwhelming^. Advancement in un- 
derstanding the fundamental physics of electron and hole 
transport in it^ and using graphene for spintronics'^ is 
critical to the realization of carbon-based electronics. 
However, to make a transition from graphene science 
to engineering and finally to a viable technology, it is 
crucial to understand the interaction of graphene with 
external parameters such as substrates^, contacts for 
measurements^, and the role of high dielectric constant 
oxides® in addition to studying the effects of intermediate 
species which are present in a particular process flow of 
graphene device fabrication. These parameters can pose 
a fundamental challenge to the ultimate realization of 
carbon electronics. In this article, we report the role of 
one of these external parameters, namely the substrates, 
in changing the electronic structure of graphene using a 
density functional based electronic structure method^'® 
and a local approximation to the exchange-correlation 
potential (LDA)2i. We use two insulating crystalline sub- 
strates, Si02 (or quartz) and AI2O3 (or sapphire), each 
with two surface terminations in our study. These are 
widely used substrates in experiments that use exfoliated 
graphene. We consider surface dangling state passivation 
with hydrogen atoms and atomic relaxations and focus 
on their role in perturbing the graphene electronic spec- 
trum. We discuss the energetics and the energy band 



structures that result from the interaction of graphene 
with these surface terminations. 

Experimentally, there are studies which attempt to un- 
derstand the atomic structure of graphene on insulat- 
ing substrateSii^, using nanometer scale microscopic tech- 
niques such as scanning tunneling and atomic force mi- 
croscopies. In addition, using Raman spectroscopy, the 
role of substrates on phonon dispersions is also reported 
which indirectly hints at the change in the electronic 
spectrum of graphene due to the presence of substrates^i. 
Moroever, we are aware of two recent density functional 
studies of monolayer graphene on a crystalline Si02 
substratei^ii^. Due to different and insufficient struc- 
tural details in those reports, we could not compare our 
results with the conclusions reached in References 12 and 
13. 

Our paper is organized as follows. In section II, we first 
discuss the computational method and the convergence 
parameters used for this study followed by the motivation 
and the procedure for building the surface models from 
their corresponding bulk counterparts. Our results of the 
effects of dangling state and its passivation on the linear 
spectrum of graphene and the role of atomic relaxation 
is discussed in section III. We describe the band struc- 
ture with the help of orbital and atom projected densities 
of states. We then present the energetics of graphene on 
both non-passivated and passivated surface terminations. 
Finally, we summarize our results and present our con- 
clusions. 
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FIG. 1: (Color online) Schematic illustration of the supercell 
structure of monolayer graphene on the top of (a) Si and 
(b) Oxygen-terminated quartz. Four unit cells of quartz are 
shown with Si atoms in blue, oxygen atoms in red and carbon 
atoms in yellow. Each Si atom is surrounded by four oxygen 
atoms. The supercell structure of graphene on the top of both 
Al and oxygen terminated sapphire is similar except that the 
successive Al and oxygen layers are stacked vertically at an 
angle. Hence we do not show them here. 



II. COMPUTATIONAL METHOD AND 
SURFACE MODELS 
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FIG. 2: Energy band structures of monolayer graphene on Si- 
terminated quartz, at high-symmetry points in the supercell 
Brillouin zone, at the interlayer distances of (a) 3 A and (b) 
2.5 A. The atoms in the supercell were not relaxed and the 
top surface was not passivated. The interlayer distance cor- 
reponds to the location of graphene plane from the topmost 
Si-plane in the Si-terminated quartz substrate. The Fermi 
energy is set at zero. The origin of occurance of linear bands 
at the F-point instead of K-point of the supercell is explained 
in the text. 



This section addresses the details of the computa- 
tional method we used followed by the procedure we 
adopted to obtain the surface models of quartz and 
sapphire from their bulk counterparts and the conver- 
gence parameters used for this study. We used a plane- 
wave based electronic structure method with local den- 
sity approximation (LDA)^ for exchange and correlation 
and the projector augmented plane-wave potential for 
electron-ion interaction-. The bulk structures of both 
the quartz and sapphire are consistent with those avail- 
able in literaturei^ and the surface models built from the 
bulk structures conform to the widely accepted a-quartz 
and Qf-alumina structures'^. We first optimized the lat- 
tice parameters of bulk crystalline quartz and sapphire. 
For both the bulk substrates, we generated hexagonal 
unit cells from the original rhombohedral cells and these 
bulk phases form layered structures (alternating cation 
and anion layers). The unit cell of quartz contains 27 
atoms with 3 Si planes and 6 oxygen planes, each plane 
containing 3 atoms, whereas the unit cell of sapphire con- 
tains 30 atoms with 4 aluminum (Al) planes and 6 oxygen 
planes again each plane containing 3 atoms. We used a 
7x7x5 k-point mesh in the hexagonal Brillouin zone 
(BZ) and a 612 eV kinetic energy cut-off. The results 
were carefully checked with respect to a larger k-point 
set and higher energy cut-offs. 



We find that the optimized lattice constants of quartz 
and sapphire are close to the experimental values (Ta- 
ble I). The in-plane and out-of-plane lattice constants 
for quartz and sapphire differ, from the experimental 
values, by less than 1 %. Using these optimized lat- 
tice parameters, we constructed surface models of both 
quartz and sapphire as follows. Four bulk unit cells were 
stacked along the c-direction (which corresponds to a 
thin-film thickness of about 22 A) and we find that 6 
X dc-c graphene, containing a total of 24 carbon (C) 
atoms (where dc-c =1-42 A) is nearly commensurate 
with the hexagonal surface of the substrates (Fig. [T|). 
The lattice mismatch of the quartz and sapphire termi- 
nations with graphene is calculated to be 0.19 % and 
0.42 % respectively. We note that the commensur ability 
of graphene with the underlying substrate is not neces- 
sary for our study since it is not focussed on epitaxial 
growth of graphene on SiC and other metal substrates. 
Therefore, the lattice mismatch values mentioned here 
serve only as an initial guideline in assessing the degree 
of distortion of strictly two-dimensional graphene in the 
presence of the substrates. In fact, we observe rippling 
of graphene layers on each of the surface terminations in 
the final relaxed structures. We note that the atomic- 
scale roughness simulated here by relaxing top few layers 
of substrates does not resemble the roughness present 
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on amorphous Si02 (or a-Si02) and AI2O3 (or a-Al203) 
surfaces that are used in experiments involving graphene. 
Simulating roughness of amorphous surfaces, using DFT 
based electronic structure method, will increase the com- 
putational burden significantly due to the requirement of 
large supercell size which is necessary to achieve rough- 
ness at the macroscopic scale. However, in the surface 
models adopted here, in their final equilibrium configu- 
rations, some roughness is present at the atomic scale. 
Periodic boundary conditions were enforced along the 
surface directions whereas a vacuum size of 10 A was 
used along the c direction to enable periodic slab calcu- 
lations. The silicon (Si) dangling states at the bottom of 
the supercell were saturated with hydrogen atoms. We 
fixed the supercell lattice parameters in all cases and only 
the atoms in the planes of the top two unit-cells of the 
substrate and the atoms in the graphene planes were al- 
lowed to relax. We used the same energy cut-off as in 
the bulk calculations but the k-point mesh in the BZ 
was chosen to be 7 x 7 x 1. For atomic relaxations, the 
total energy was assumed to have converged when all the 
components of the Hellman-Feynman forces were smaller 
than 0.01 eV/A. 
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FIG. 3: Energy band structures of monolayer graphene on 
Al-terminated sapphire, along high-symmetry points in su- 
percell Brillouin zone, at interlayer distances of (a) 4.0 A and 
(b) 2.7 A. The atoms in the supercell were not relaxed and 
the top surface was not passivated. The interlayer distance 
correponds to the location of the graphene plane from the 
topmost Al-plane in Al-terminated sapphire substrate. The 
Fermi energy is set at zero. 



III. ENERGY BAND DISPERSIONS, 
DENSITIES OF STATES AND ENERGETICS 

In this section, we discuss the role of dangling states, 
their passivation and the role of atomic relaxations on 
the electronic band structure of graphene on the sur- 
face terminations we considered here. We note that both 
the surface terminations of quartz possess large number 
of dangling states whereas in sapphire, by construction 
of the surface models from bulk, the number of nearest 
neighbors of both Al and oxygen is found to be optimal. 
Therefore, for passivation studies, only quartz surfaces 
were considered. We passivated the dangling states with 
hydrogen atoms. There are four different possibilities 
that emerge from considering whether surface dangling 
states are passivated or not and whether or not atomic 
relaxations are performed. 

We first discuss Si-terminated quartz. We did not pas- 
sivate the dangling states and did not relax the atomic 
positions. To get a rough estimate of the interlayer sep- 
aration at which graphene electronic structure is mini- 
mally perturbed, we placed the graphene layer at vari- 
ous distances from the top substrate layer. The choice 
of these distances and distances chosen for other surface 
terminations, are arbitrary. We chose bilayer graphene 
interlayer distance (which is 3.34 A) as a reference. At 
a distance of 3 A and above, from the top Si surface, 
graphene retains the linear spectrum but at a smaller 
distance (of 2.5 A), the linear bands are perturbed (Fig. 
[2l[a)and (b)). The occurance of the linear spectrum at 
the F-point of the supercell BZ instead of the K-point is 
due to crystal symmetry of the supercell and its relation 
to the monolayer graphene lattice symmetry. The origin 
of such shifting of the location of the electron-hole sym- 
metry from the K-point to the F-point is seen in calcula- 
tions involving sub-nionolayer alkali metal adsorption on 
graphene surfaces^i. We discuss the origin of this shift 
in the Appendix section. On relaxing the atomic posi- 
tions of the top few layers of the quartz substrate and 
graphene, we get an equlibrium distance of 3 A between 
the top substrate layer and graphene due to which the 
perturbation to electronic spectrum is minimal, as in the 
non-relaxed case (Figure not shown). 

We now discuss Al-terminated sapphire and its inter- 
action with graphene. For the non-passivated surface 
without atomic relaxations, the dispersion curve shows 
the dangling states at the Fermi level are resonance with 
the carbon pz orbital at interlayer distances as high as 
4 A and as low as 2.7 A (Fig. [3)Ja) and (b)). However, 
with atomic relaxations, we get an equlibrium separation 
of 2.7 A and a similar energy dispersion curves (Figures 
now shown). 

The situation is similar in both oxygen-terminated 
quartz and sapphire. On non-passivated oxygen- 
terminated quartz without atomic relaxations, at dis- 
tances of 1.76 A and 2.5 A, graphene features appear in 
the band structure (Fig. Ufa) and (b)) but the electron- 
hole symmetry is located above the Fermi level but an 
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FIG. 4: Energy band structures of monolayer graphene on 
oxygen-terminated quartz, at a interlayer distance of (a) 1.76 
A and (b) 2.5 A. The atoms in the supercell were not relaxed 
and the top surface was not passivated. The interlayer dis- 
tance correponds to the location of the graphene plane from 
the topmost oxygen plane in oxygen-terminated quartz sub- 
strate. The Fermi energy is set at zero. 



energy gap opens up for the case of 2.5 A. Opening of 
an energy gap may be due to the breaking of sub-lattice 
symmetry in graphene as observed in epitaxial growth of 
graphene on SiC substrate— but is most likely due to ap- 
proximations in DFT. On relaxing the atomic positions, 
oxygen-terminated surfaces in both quartz and sapphire, 
maintain an equilibrirm separation which is quite close 
to the graphene layer. As a result, we find a strong per- 
turbation to the linear spectrum of graphene (Fig. 
and (b)). 

When we add a second carbon layer on oxygen- 
terminated quartz and relax the atomic positions, we re- 
cover linear spectrum of graphene (Fig. EJc)). This hints 
at the existence of a buffer layer in studies of graphene 
on oxygen-terminated quartz. However, our calculations 
suggest more than two carbon layers are needed to re- 
cover graphene linear energy spectrum in case of Al- and 
oxygen-terminated sapphire. 

These equlibrium distances, for each non-passivated 
surface terminations with atomic relaxations, are listed 
in Table I. It also lists interatomic distances between the 
Carbon and the corresponding atom of the surface ter- 
minations in the bulk or molecular phase. Since all the 




FIG. 5: Energy band structures of monolayer graphene on (a) 
oxygen-terminated quartz, (b) oxygen-terminated sapphire, 
along high-symmetry points in supercell Brillouin zone, com- 
puted at the respective equilibrium separations shown in Ta- 
ble I and (c) two carbon layers on the top of the oxygen- 
terminated quartz seems to restore the linear behavior at the 
Fermi level. The atoms in the supercell were relaxed but the 
top surface was not passivated. The Fermi energy is set at 
zero. 



equilibrium distances are larger than that in the corre- 
sponding interatomic distances in the bulk or molecular 
phases, it suggests that graphene does not form a bulk- or 
molecular-like phase with the surface terminations used 
in this study. 

We now discuss the effect of passivation of dangling 
states on the band structure of graphene. Since graphene 
on non-passivated Si-terminated quartz already shows 
the a linear spectrum, we considered passivating only the 
dangling states of oxygen-terminated quartz. We did not 
relax the atomic positions. Linear bands seem to ap- 
pear above the Fermi level for distances below 2 A but 
right at the Fermi level in case of distances above 2 A 
(Fig. IHKa) and (b)). It is likely that relaxing the atomic 
positions may retain the graphene features but with a 
different equilibrium separation than the non-passivated 
surface. We did not consider this case explicitly in our 
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calculations. 

Our results suggest that presence or absence of dan- 
gling states is not as effective in distorting the linear 
band structure as the atomic relaxations. It is the equi- 
librium distance that dictates the survival of graphene 
features in the dispersion spectrum. However, a case in 
which relaxations of passivated substrate and graphene 
provide a favorable equilbrium distance for retention of 
the graphene features, is not ruled out. 

To understand the perturbations to the linear spec- 
trum of graphene in case of non-passivated relaxed sur- 
face, we plot atom and orbital projected densities of 
states (or DOS) (Fig. The top panel (a) shows DOS 
for oxygen terminated quartz, projected on to s- and p- 
states of oxygen atom of the top surface layer and carbon. 
We see a strong hybridization between oxygen-p and C-p 
orbitals in the vicinity of the Fermi level. A similar con- 
clusion was reached in a recent DFT-based calculation 
on monolayer graphene on oxygen-terminated quartzi^. 
Such resonance structures are also seen in the DOS of 
both Al- and oxygen-terminated sapphire (panels (b) and 
(c)). This explains why graphene 7r-orbitals, after atomic 
relaxations, cannot retain their linear dispersion in the 
presence of non-passivated surfaces. We also estimated 
average height fluctuations of two-dimensional graphene 
due to the presence of substrate and the atomic relax- 
ations of the top surface layers. In all cases, we find that 
the deviations are ~ 0.05 A which is small compared to 
suspended graphene^^. 

Table I shows the binding energy values for graphene 
on the four surface terminations considered in this study. 
In the most relevant case i.e. no surface passivation but 
with atomic relaxations, these values hint at non-bonding 
nature of graphene to the underlying oxide substrates. 
The binding energy values are obtained by using the fol- 
lowing definition. 



^bind ^ ^(supercell) - .B(Gr) - £;(substrate) (1) 

where i?(supercell) denotes the total energy of the su- 
percell containing the substrate and a graphene layer. 
E{Gr) and i?(substrate) denote, respectively, the total 
energies of isolated graphene and isolated substrate in 
the same supercell set-up, with the same energy cut-off 
and k-point mesh as that of the combined graphene and 
substrate calculations. 

However, we find that passivated but non-relaxed sur- 
faces are energetically more favorable compared to both 
non-passivated non-relaxed as well as non-passivated 
relaxed surfaces. This suggests an existence of a 
metastable configuration of graphene in presence of the 
non-passivated surfaces (both with and without atomic 
relaxations). Of course, passivated and relaxed surfaces, 
which we did not explicitly calculate here, can be lower 
in energy than any of these cases. 



TABLE I: Lattice parameters (in A) of bulk quartz and sap- 
phire, average interplanar distances (in A) of graphene from 
the four underlying surface terminations and their binding 
energies (in eV/atom), in case of non-passivated but relaxed 
surfaces. The numbers in parenthesis are the out-of-plane lat- 
tice constants of the bulk phases and those in square paren- 
thesis are interatomic distances between Carbon and the cor- 
responding atom of the surface terminations in the bulk or 
molecular phase. It should be noted that, for isolated sub- 
strate calculations, we passivated the top layer with hydrogen 
atoms to keep the supercell non-magnetic. 

Lattice Paramters d(C-x) E""'' 

(x = Si,Al,0) 

Bulk quartz 

This work 4.914 (5.408) 

Expt." 4.913 (5.405) 

Si-terminated 3.0 (1.89)' 10.061 

Oxygen-terminated 1.76(1.3)'= 0.581 



Bulk sapphire 

This work 4.907 (4.908) 

Expt."* 4.943 (4.907) 

Al-terminated 2.7 (1.89-2.19)" 1.017 

Oxygen-terminated 2.15 0.689 



"Reference 14, 16 
''Reference 17 
"^Reference 18 
''Reference 15, 19 
"^Reference 20 



IV. SUMMARY AND CONCLUSION 



In summary, we have studied, using a first princi- 
ples DFT method, the effect of two crystalline insulat- 
ing substrates, quartz and sapphire, on the electronic 
structure of monolayer graphene. We considered the ef- 
fect of surface passivation and atomic relaxations on the 
linear spectrum of graphene. Dangling states seems to 
be ineffective in distorting the linear band structure of 
graphene except shifting the spectrum above Fermi level. 
It is the atomic relxations which dictate the equlibrium 
separation between graphene and the topmost surface 
layer and this distance decides the strength of perturba- 
tion on linear bands. Si-terminated a-quartz retains the 
graphene band structure on non-passivated surface even 
with atomic relaxations. This leaves oxygen-terminated 
quartz surface where atomic relaxations play an impor- 
tant role. Both Al- and oxygen-terminated sapphire per- 
turb the linear bands whether or not their dangling states 
are passivated and whether or not atomic relaxations 
are performed. Two carbon layers are necessary to re- 
cover the linear band structure of graphene on the top of 
oxygen-terminated quartz whereas our calculations hint 
at more than two carbon layers are needed in case of Al- 
and oxyen-terminated Sapphire. Energetically, graphene 
on non-passivating surface terminations is metastable. 
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FIG. 6: (Color online) Same as Figure 4 but now for passi- 
vated surface of oxygen terminated quartz. 
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Appendix 

Quartz and sapphire substrates and the graphene 
layer, all have a hexagonal unit cell. Let Ai, A2 and 
A3 be the lattice vectors of the entire graphene and sub- 
strate system and Bi, B2, and B3 be the correponding 
reciprocal lattice vectors. Similarily, let £Z2 , lis) and 

{ bi, 62 and 63) be the triad of primitive and reciprocal 
lattice vectors of the graphene layer alone. 

We note that the lattice structure of the graphene plus 
substrate system is a (2 y^S) x 2 ■y/(3))R30° recon- 
struction of the lattice structure of the graphene layer 
alone. Thus A = 2y/{3)a where A is the lattice con- 
stant of the entire system and a is the lattice constant 
of the two-dimensional graphene layer. Since both unit 
cells are hexagonal this implies that the reciprocal lat- 
tice structure of the graphene layer on top is a {2^/{3) x 
2-\/(3))R30° reconstruction of the reciprocal lattice struc- 
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FIG. 7: (Color online) Atom and orbital projected density 
of states for a C atom in monolayer graphene and the atoms 
of the topmost planes in (a) oxygen-terminated quartz (b) 
Al-terminated sapphire and (c) oxygen-terminated sapphire 
substrates computed at the equilibrium separations shown in 
Table I. The atoms in the supercell were relaxed but the top 
surface was not passivated. The Fermi energy is set at zero. 



ture of the graphene plus substrate system below. Taking 
^1 to be along x, 



oi = ax 



a2 = o(-x+ — ^y) 



li=2y(3)a(^x+iy) 
= 2 7(3) ay 



(A.l) 
(A.2) 

(A.3) 
(A.4) 
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As = cz 



(A.5) 



t 1 

oi = — (x 7=— y) 

a 7(3) 



(A.6) 



(A.7) 
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FIG. 8: (Color online) Schematics of supercell Brillouin zone 
(small) superposed on the graphene honeycomb Brillouin zone 
(large). The reciprocal lattice vectors are shown along with 
the overlap of F to K vector of graphene with the Bi vector 
of the supercell. The overlap is due to the symmetry reasons 
and is explained in the text. 

The reciprocal lattice vectors are, 



Bo = 



27r 



2^(3) a 7(3) 



27r 



(- 



2^(3) a 7(3) 



-y) 



(A., 



(A.9) 



Figure 8 shows the BZ set-up with these vectors using 
equations A.6 to A.9. 

Now if we take the same F point as the center of both 
the Brillouin zones and draw the reciprocal unit cell of 
both the structures we realize that the K point of only 
the graphene layer lies on the reciprocal lattice vector 
_Bi of the graphene plus substrate system (as shown in 
figure). This means that the K point of graphene layer 
folds in by symmetry onto the F point of the entire sys- 
tem. 
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